Collisional properties of trapped cold chromium atoms 
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We report on calculations of the elastic cross section and thermalization rate for collision between 
two maximally spin-polarized chromium atoms in the cold and ultracold regimes, relevant to buffer- 
gas and magneto-optical cooling of chromium atoms. We calculate ab initio potential energy curves 
for Cr2 and the van der Waals coefficient Ce, and construct interaction potentials between two 
colliding Cr atoms. We explore the effect of shape resonances on elastic cross section, and find that 
they dramatically affect the thermalization rate. Our calculated value for the s-wave scattering 
length is compared in magnitude with a recent measurement at ultracold temperatures. 



Collisions of atoms at ultracold temperatures have re- 
ceived considerable attention because of their importance 
in cooling and trapping of atoms and their role in high 
precision spectroscopy and Bose-Einstein condensa- 
tion Q. Processes occurring in the cold regime are sen- 
sitive to the details of the interaction potentials between 
the colliding systems over an extended range of inter- 
nuclear separations. Recent experiments with chromium 
^ IE IE Q emphasize the need for theoretical studies of 
Cr scattering properties. The interest in cooling Cr stems 
from its particular properties; in its ground state ^^3, it 
possesses a very large magnetic moment, 6 /is (/is^ Bohr 
magneton) , making it an ideal atom for buffer-cooling in 
a purely magnetic trap , as well as for magneto-optical 
trapping P|. In addition, anisotropic long-range interac- 
tions, such as chromium's magnetic dipole-dipole inter- 
actions, may lead to novel phenomena in BECs 0,^3- 
The existence of a stable fermionic isotope, '^'^Cr, opens 
the possibility of obtaining fermionic degenerate gas us- 
ing sympathetic cooling. Chromium was also used in a 
new cooling scheme Q , where ultracold chromium atoms 
were loaded from a MOT into an loffe-Pritchard mag- 
netic trap, and cooled below IQQjjLK. However, a not 
so-desirable byproduct of large-spin collision is inelastic 
"bad" scattering rates that deplete the trap. 

On the theoretical front, the electronic spectrum of 
the Cr2 dimer poses a considerable numerical challenge. 
Chromium is the first atom in the periodic table with a 
half-filled d shell (the ground electronic configuration is 
Cr(3(i^4s,'' 5*)) and the Cr2 dimer is one of the most ex- 
treme cases of multiple metal-metal bonding. To date, 
the best attempt to calculate its interaction potential 
curves is a multiconfiguration second-order perturbation 
theory with complete active space self-consistent field 
(CASSCF/CASPT2) [ll|l3. While some information 
on the spectroscopy of the ground electronic state 
molecular symmetry) exists, there is no data available for 
the interaction of two maximally spin-stretched Cr atoms 
in the ^^S^ molecular symmetry, i.e. total spin, 5 = 6. 

In this communication, we explore the collisional prop- 
erties of Cr atoms in the cold and ultracold tempera- 
ture regimes by revisiting the electronic structure of the 
dimer. More accurate Born-Oppenheimer potential en- 



ergy curves dissociating to two ground Cr atoms were 
obtained. The van der Waals interaction coefficient, 
Cg, was obtained semi-empirically from available bound 
atomic transition matrix elements and photoionization 
cross sections to be Ce = 745 ± 55 a.u., where 1 a.u. — 
9.57 X 10~*° Jm^. The elastic cross section and colli- 
sion rate coefficient were calculated using the newly con- 
structed potential curves and compared against two re- 
cent measurements of elastic rate coefficients '5','?|. We 
also investigated the bound and resonance structure in 
the ^"^5]+ potential energy curve and determined the res- 
onance positions and widths as a function of rotational 
angular momenta. 

The potential curves for Cr2, shown in Fig. ^ were 
constructed from three regions joined smoothly together. 
We first computed ab initio potential curves using the 
CASSCF/CASPT2 method III El- The CASSCF wave 
function is formed by distributing 12 electrons in the 
3d and 4s derived active orbitals while keeping the in- 
active Is, 2s, 2p, 3s, and 3p derived orbitals occupied. 
The remaining dynamical electron correlation energy 
is obtained through second order perturbation theory 
(CASPT2). The basis set used in the calculations is of 
the atomic natural orbital (ANO-RCC) type contracted 
to 9s8p7d5f3g. This basis set is relativistic and includes 
functions for correlating the 3s and 3p electrons |3| • The 
Douglas-KroU Hamiltonian was used with Fock-type cor- 
rection 0.5 * gi, see |lj|. The full counterpoise method 
was used to correct energies for the basis set superpo- 
sition error (BSSE). Convergence to 10^"'^°, in hartrees, 
was achieved, and numerical accuracy in computed bind- 
ing energies was about 10^^. 

For separations R < R^, where R^ is the smallest sepa- 
ration of the ab initio data for the potential energy Vx (i?) , 
each curve was joined smoothly to the exponential form 
V\{R) — c\exp{—bxR), with the coefficients cx and bx 
determined by matching both the potential curve and its 
first derivative continuously at i?^. At large values of R, 
the ab initio data were matched to the asymptotic form 

Vx{R)^-^+AxR-'e'^" , 
where Cq is the van der Waals coefficient and the 
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FIG. 1: Potential curves for the ground state manifold of Cr2 
computed with the CASSCF/CASPT2 method. The max- 
imally spin-streched electronic state, ^"^S^, is shown in red 
and in detail in the inset. This state supports 30 vibrational 
levels. 



parameters of the exchange energy are determined ac- 



gyar 

cording to Smirnov and Chibisov |l6j : = ^ ~ 1 &iid 
(3 = 21, where I is the ionization energy of the atom 
(/ = 0.248664314 a.u. for Cr). The parameters Ax were 
found by fitting the ab initio curves at separations where 
the exchange energy was still considerable (e.g. R be- 



tween 10 and 14 a.u. 
calculated using 



for 
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The Cfi coefficient was 




where contributions of bound-bound, bound-free, and 
free- free transitions are given by the first, second, and 
third term, respectively. This expression is derived from 
the London's [l7| formula, assuming Vi = 1 — Ei /Eq and 
e ~ ~E/Eq, where Eq, Ei, and E are the ground, the i^^ 
excited state and continuum energies, respectively, foi 
are the oscillator strengths pertaining to transitions to 
the ground states, and dfoE accounts for the continuous 
spectrum. The auxiliary function G{z) is given by 



G{z) 



{df/de)de 
(l + e)(2 + e) 



(2) 



Values of oscillator strengths and energy levels for dis- 
crete transitions were taken from the NIST Atomic Spec- 
tra Database and Ref. . The continuous oscillator 
strength, df/de, was found using Verner's ^| analytic 
fits for partial photoionization cross sections. If we mea- 
sure dE in atomic units then 
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(3) 



where a is the fine-structure constant, ao is the Bohr 
radius, and CTph is the photonization cross section. The 
dimensionless df/de is obtained by multiplying Eq.ljSJ by 
\Eo\ (also in a.u.). 

We assess the quality of the semi-empirical computa- 
tion of the van der Waals coefficient by satisfying two 
sum rules, namely, the zero-th and the inverse second 
moments. 



^(-2) = ao 




(4) 



where N is the number of electrons (24 for Cr), and ao is 
the static polarizability. We obtained J^i foi — 22.3 and 
ao — 85.00 a.u., in agreement with a recommended value 
82±2G%a.u. |23|). The resulting dispersion coefficient is 
expected to have an accuracy of about 7%, Cg — 745 ±55 
a.u. 

Using these potential curves, we computed the elastic 
cr^j cross section for the maximally spin-aligned molecu- 
lar state, i. e. the state. [Results for other molec- 
ular states and also inelastic and spin relaxation cross 
sections will be provided in a future publication.] The 
elastic cross section for the collision of two ^^Cr atoms, 
composite bosons, expanded over the rotational quantum 
number, /, is [23 



8n 
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^(2Z + l)sin2 5i^ 



(5) 



where E = H k is the kinetic energy of relative mo- 
tion, ^ is the reduced mass, and 5i{k) is the l-ih scatter- 
ing phase shift in electronic state A. In the low-energy 
limit, the elastic cross section behaves as 
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where scattering length ax is determined at the zero- 
energy s-wave limit, ax = — lim^^o tan i5g (/c) , while 
the effective range rx can be found by fitting the cross 
section, from an i nteg ral expression |2^|. or using quan- 
tum defect theory pa|. 

The cross section as a function of the collision energy 
is shown in Fig.|21for three different interaction po- 

tentials; each constructed by matching at large distances 
to the upper, lower, and mean values of Cg. Note that 
at collision energies larger than E ^ 10~^ a.u., aside 
from the shape resonance structure, the cross sections 
differ little in magnitude, but are dramatically different 
in the s-wave limit. We find the scattering length for 
the state, 013, to be large and negative (see Ta- 

ble P), indicating that evaporative cooling should be efR- 
cient, but that a Bose-Einstein condensate of ^^Cr would 
not be stable for large number of atoms. This contrasts 
Ref. jsl] in which the extracted value for the scattering 
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FIG. 2: Elastic cross sections for collision between two 
bosonic Cr atoms as a function of energy in the cold and 
ultracold regimes. The three curves are labeled by the value 
of the Ce coefficient used to construct the interaction poten- 
tial in each case. See text for details. The influence of shape 
resonances is dramatic. 



length, obtained from a fit to the experimental collision 
rate coefficient, Fig. 13 agrees in magnitude with our up- 
per bound calculated value, but differs in the prediction 
of the sign (see also the discussion of the rates below). 
Our calculated effective range expansion coefficients also 
agree very closely with the results obtained using the 
quantum-defect method of Gao 23] . 

As the collision energy increases, the appearance of 
shape resonances can lead to enormously large cross sec- 
tions. Althougth the details of the resonance structure, 
i. e. profile and energy position, depend on the details 
of the potentials, the overall effect on the coUisional rate 
coefficient is only minimal at temperatures for which the 
shape resonances matter. This is portrayed in Fig. |2| 
for the ultracold and cold temperatures. The rate coef- 
ficients for the three different potentials are similar for 
temperatures higher than 100 mK. Our rate coefficients 
are larger than those measured by J] by more than an 
order of magnitude. We do not yet understand the origin 
of this discrepancy- an attempt at modifying the poten- 
tial in the long-range and in the short-range resulted in 
practically similar results. The shape and magnitude of 
the calculated rate coefficient for Ce = 800 a.u. agrees 
well in both magnitude and shape with the MOT exper- 
imental results, see Fig. 3 of Ref. The decline in 
the rate coefficient for T > 100 /iK, as seen in our cal- 
culation, followed by a rise for T > 1 mK, is due to the 
appearance of the I = 4 shape resonance in Fig. |2| The 
experimental determination of the dip in the elastic rate 
could be used to infer a more accurate value of Cg and 
hence the scattering length as well. 

In Fig. ^ we give the bound and shape resonance 
structure in the ^'^S^ potential energy curve with Cg = 
745 a.u., as a function of the rotational angular momen- 
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FIG. 3: The calculated collision rate coefficient as a function 
of temperature. At temperatures above 100 mK, the different 
calculated rates are practically the same. The plateau region 
is shown in the inset. 

tum, I. Due to symmetry, only even values of I appear 
in the collision. There are 30 bound vibrational levels 
for the I ~ partial wave. The grouping of the levels 
in solid and dashed lines indicates the appearance of an 
additional shape resonance. Table UTI gives our calculated 
position and width for a number of rotational shape reso- 
nances. The narrowest resonance is for the I = 19 partial 
wave with a lifetime of about two seconds. This reso- 
nance, however, is not populated in a bosonic collision. 

A plateau in the observed elastic collision rate coeffi- 
cient in the range 10 mK < T < 300 mK exists 0, Q 
and is reproduced in our calculation in Fig. 13 It appears 
that this plateau is produced by the confluence of many 
coUisionally-excited shape resonances in this temperature 
range. 

We have calculated ab initio interaction potentials for 
collision of two Cr atoms, and obtained the long-range 
van der Waals coefficient for the first time. Elastic colli- 
sion cross section and rate coefficient for two maximally- 
stretched chromium atoms have been computed in the 
cold and ultracold regimes and compared with available 
experiments. The effect of partial-wave shape resonances 
is studied. The s-wave scattering length and effective 
range have been obtained in accord with experiment. 
The disagreement with the buffer-gas cooling rate coeffi- 
cient results is not fully understood and the inelastic loss 
rate will be studied in the near future. The main loss 
rate mechanism is expected to be spin dipole interaction. 
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FIG. 4: The bound and shape resonance structure in the Cr2 

^^E+ potential energy curve as a function of the rotational TABLE I: The calculated scattering length and effective 
quantum number I. range for different interaction potentials, all in atomic units. 
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TABLE II: The rotational shape resonance parameters, po- 
sition and width, in the ^^'S^ potential. Note that only even 
partial waves contribute to the scattering of spin polarized 
bosons (odd partial waves are included for completeness). 



